function [upper,lower,Insigma] = theta1BoundsInsigma(nm,nn)

%% Uper bounds
upper.alpha = 10; 
upper.rhor  = 1;
for i=1:2*nm+nn
   for j=1:2*nm+nn
        name = ['phi',num2str(i),num2str(j)];
        upper.(name) = 10;
   end
end
for i=1:2*nm+nn
    name = ['h0Q',num2str(i)];
    upper.(name) = 10;
end
for i=1:2*nm+nn
   for j=1:i
        name = ['sigma',num2str(i),num2str(j)];
        if i == j
            upper.(name) = 0.05*10;
        else
            upper.(name) = 0.5;
        end
   end
end

%% Lower bounds
lower.alpha = -10; 
lower.rhor  = 0.0;
for i=1:2*nm+nn
   for j=1:2*nm+nn
        name = ['phi',num2str(i),num2str(j)];
        if i == j
            lower.(name) = 1D-6;
        else
            lower.(name) = -10;
        end
   end
end
for i=1:2*nm+nn
    name = ['h0Q',num2str(i)];
    lower.(name) = -10;
end
for i=1:2*nm+nn
   for j=1:i
        name = ['sigma',num2str(i),num2str(j)];
        if i == j
            lower.(name) = 1D-6;
        else
            lower.(name) = -1;
        end
   end
end
%% Insigma
Insigma.alpha = 0.01; 
Insigma.rhor  = 0.1;
for i=1:2*nm+nn
   for j=1:2*nm+nn
        name = ['phi',num2str(i),num2str(j)];
        Insigma.(name) = 0.01;
   end
end
for i=1:2*nm+nn
    name = ['h0Q',num2str(i)];
    Insigma.(name) = 0.05;
end
for i=1:2*nm+nn
   for j=1:i
        name = ['sigma',num2str(i),num2str(j)];
        Insigma.(name) = 0.01;
   end
end

end

